Solving mazes with memristors: a massively-parallel approach 
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Solving mazes is not just a fun pastime: they are prototype models in several areas of science and 
technology. However, when maze complexity increases their solution becomes cumbersome and very 
time consuming. Here, we show that a network of memristors - resistors with memory - can solve 
such a non-trivial problem quite easily. In particular, maze solving by the network of memristors 
occurs in a massively parallel fashion since all memristors in the network participate simultaneously 
in the calculation. The result of the calculation is then recorded into the memristors' states, and 
can be used and/or recovered at a later time. Furthermore, the network of memristors finds all 
possible solutions in multiple-solution mazes, and sorts out the solution paths according to their 
length. Our results demonstrate not only the first application of memristive networks to the field of 
massively-parallel computing, but also a novel algorithm to solve mazes which could find applications 
in different fields. 

PACS numbers: 87.18. Sn, 73.63.-b, 73.50.Fq 



I. INTRODUCTION 

Mazes are a class of graphical puzzles in which, given 
an entrance point, one has to find the exit via an intri- 
cate succession of paths, with the majority leading to a 
dead end, and only one, or few, correctly "solving" the 
puzzle. Mazes - sometimes also called labyrinths - have 
been known since ancient times, the oldest presumably 
being the one created by Daedalus in Crete, as passed 
on by Greek mythology a few thousand years ago. They 
are used as prototype models in graph theory, topology, 
robotics, traffic optimization, psychology, and in many 
other areas of science and technology Uj-^. The ability 
to solve a maze is particularly important for transporta- 
tion and robot control systems, e.g., to find the shortest 
path from a given point to another. In neuroscience, dif- 
ferent maze solving algorithms could be compared with 
the approaches used by humans to solve mazes. In this 
way we could better understand the functioning of the 
human brain and advance new modeling of neural pro- 
cesses. 

Algorithms to solve mazes vary from the simplest - 
and extremely slow - "random mouse" to mathematical 
search algorithms that operate on a sequential fashion 
to find the exit. However, all these methods suffer from 
very slow solution times when the complexity of the maze 
increases, with solution times sometimes increasing dra- 
matically with increasing local connectivity of the net- 
work. Moreover, some researchers have demonstrated 
that certain biological and chemical systems can solve 
mazes [H |9]. For instance, in Ref. [4J an unexpected 
behavior of a primitive organism was revealed by show- 
ing that an amoeba finds the minimum-length path be- 
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tween two points in a labyrinth connecting separate food 
sources. Such methods, however, are also slow and do 
not seem to be suitable for large mazes with complex 
connectivity. 

In this paper, we suggest and demonstrate a new strat- 
egy for solving mazes that is instead based on massively- 
parallel computation as afforded by a network of mem- 
ristors (short for memory resistors) [lOj. Memristors are 
resistors whose resistance depends on the state history of 
the system, and can therefore record their past dynam- 
ics. These systems, which belong to the larger class of 
memory circuit elements - that includes also memcapaci- 
tors and meminductors [llj - are attracting considerable 
attention due to their usefulness in diverse research ar- 
eas [12j ranging from memories p!3l [T4] to neuromorphic 
computing and learning p!5UT7) . 

A network of memristors, complemented by some other 
standard electronic elements (such as field effect transis- 
tors), forms a memristor processor that, at the present 
level of technology, can be fabricated experimentally. In 
this work, instead, we use numerical simulations to ex- 
plicitly show that a memristor processor is able to solve 
mazes. Most importantly, unlike existing approaches, the 
memristor processor requires only one step to find the 
maze solution. We emphasize though that the type of 
massive parallelism used in the memristive processor can 
be thought of as analog parallelism. It is essentially differ- 
ent from that used in conventional computers, and rather 
shows some similarity with massively parallel computing 
with organic layers demonstrated recently [18j. 

It is also worth noting here that we could perform maze 
solving with a network of meminductors or memcapaci- 
tors. We choose to work with memristors since they are 
the most studied so far, and can be easily realized exper- 
imentally, thus allowing a practical and straightforward 
implementation of our algorithm. 

The remainder of this paper is organized as follows. 
In Sec. [TTj we introduce the main part of our proposal. 



FIG. 1: (Color online) Maze mapping into a network of memristors (a.k.a. memristive processor). Right panel. The maze is 
covered by an array of vertical and horizontal lines having the periodicity of the maze. Left panel Architecture of the network 
of memristors in which each crossing between vertical and horizontal lines in the array (in the right panel) is represented by a 
grid point to which several basic units consisting of memristors and switches (field-effect transistors) are linked in series. The 
maze topology is encoded into the state of the switches such that if the short line segment connecting neighboring crossing 
points in the array crosses the maze wall then the state of the corresponding switch is "not connected" (shown in red). All 
other switches are in the "connected" state. The external voltage (V) is applied across the connection points corresponding to 
the entrance {V) and exit (ground, GND) points of the maze. 



namely, the memristive processor. We discuss all impor- 
tant operation details of the memristive processor includ- 
ing its initialization, maze mapping and, finally, finding 



the maze solution. Sec. Ill illustrates dynamics of maze 
solving by the memristive processor. In this section, we 
first introduce a mathematical model of the unit part of 
the processor, the memristor. Numerical simulations of 
solutions of both single-path and multiple-path mazes by 
the memristive processor are presented. It is important 
to note that in the case of multiple-path mazes, the mem- 
ristive processor not only finds all solutions but also sorts 
them out according to their length. We conclude with 
some final remarks and considerations on the proposed 
approach. 



II. MEMRISTIVE PROCESSOR 

Let us then start by mapping a given maze into a net- 
work of memristors. This is shown in Figure [1] First 
of all, we superimpose periodical arrays of vertical and 
horizontal lines on the maze. The period of this ar- 
ray corresponds to the intrinsic period of the maze (for 
non-periodic mazes, the period of line arrays should be 
selected in such a way to take into account all impor- 
tant maze features). The crossing points of vertical and 
horizontal lines define grid points of the memristive net- 
work. The network consists of basic units (memristors 
plus switches) connecting grid points. Since the direction 



of current flow in the network is not known a priori, the 
polarity of adjacent memristors (indicated by the black 
thick line in the memristor symbol in Figure T]) is chosen 
to be alternating. It is assumed that external signals can 
be applied to any grid points for the purpose of initial- 
izing the memristors' states as well as to read the cal- 
culation results. The externally-controlled switches are 
used to define the topology of the maze: a maze wall is 
modeled by a switch in the "not-connected" state. Such 
an architecture allows modeling different mazes on the 
same memristive processor without the need to fabricate 
a specific processor for each maze. 

The processor initialization can be done by simultane- 
ous application of GND, V\ voltages in a chessboard-like 
pattern to all grid points of the memristive network (see 
Fig. |2|i). The calculation performed by this memristive 
processor occurs when a constant voltage V is applied 
across the two grid points corresponding to the entrance 
and exit points of the maze as shown in Figure [T] In this 
case, the current flows only along those memristors that 
connect the entrance and exit points. The state of these 
memristors is changed by the current, thus the maze is 
solved in a massively parallel way, since all memristors 
in the network participate simultaneously in the calcula- 
tion. Specifically, assuming that by the initial moment of 
time all memristors were initialized in the high resistance 
("OFF") state, as time passes, every other memristor 
along the solution path changes its resistance, eventually 
switching into the low-resistance ("ON") state. There- 
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FIG. 2: (Color online) a) Initialization of the network of mem- 
ristors can be performed by a simultaneous application of 
GND and appropriately selected Vi voltages in a chessboard- 
like pattern to all grid points for a sufficiently long period 
of time to securely switch memristors to the "OFF" state. 
Note that since the memristors' polarities are alternated in 
the network, the corresponding voltage polarities are also al- 
ternated. During the initialization, all switches should be 
in the "connected" state or should already encode the maze 
(as shown on the plot), b) In order to read a memristor's 
state, such a memristor (shown in the center) can be discon- 
nected from the rest of the circuit putting the corresponding 
switches to the "not-connected" state (shown in read) and 
its own switch to the "connected" state and tested by the 
application of a short small-amplitude voltage pulse (or dou- 
ble bi-polar voltage pulse to minimize the disturbance of the 
memristor's state) of an amplitude V2. 



fore, the chain of memristors in the "ON" state (or in 
an intermediate state if memristors did not have enough 
time to completely switch to the "ON" state) connecting 
the entrance and exit points represents the maze's solu- 
tion. A possible approach to read the calculation result 
is described in Fig. ^p. 

Furthermore, the states of the memristors connecting 
the entrance and exit points represent a solution at any 
given finite time after the initial one. If there are multiple 
paths, then the shortest one would contain less memris- 
tors and thus offers less resistance than the longest one, 
with all intermediate paths (in terms of length) offer- 
ing a proportionate resistance. Therefore, since current 
flows in inverse proportion to the resistance of a path, 
at any given time, the change of state of a given path 
is proportional to the current in the path. The different 
paths of the maze can then be identified by the different 
state their memristors have during (or after) the switch- 
ing process. In the next section, we explicitly discuss this 
solution feature. Here, we only point out that the mem- 
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FIG. 3: (Color online) Solution of a single-path maze, a. 
Network state at t = 0.075s. The chain of memristors in the 
low-memristance state (shown by red dots) clearly connects 
the entrance and exit points of the maze (note that memris- 
tors in the "OFF" state - when Rf^^ (t = 0) > 90 Ohms - are 
not shown). Here, every other memristor along the solution 
path is in the low-memristance state, b. Network state at 
t — 0.12s. Note that at t = 0.1s the sign of the applied volt- 
age has been changed. At t = 0.12s, each memristor along the 
solution path shows the maze solution. The resistance is in 
Ohms, the voltage is in Volts, and the current is in Amperes. 



ristive processor can be reduced to a network of standard 
resistors and, in principle, this could also solve a maze. 
However, in this case, additional external memory and lo- 
cal voltage measurement hardware will be required as the 
resistive network does not store the calculation results. 
In addition, in the case of a standard resistive network, 
previous calculation results can not be easily used in sub- 
sequent calculations. Finally, a network of memristors is 
a complex non-linear dynamical system that can be po- 
tentially used to solve a larger class of problems, of which 
optimization ones - as those described in this paper - are 
only a small class. 

Experimentally, the suggested network could be fab- 
ricated using, e.g., CMOL (Cmos+MOLecular-scale de- 
vices) architecture [19j combining a single memris- 
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tor layer with a conventional CMOS (complementary 
metal-oxide-semiconductor) layer. It was demonstrated 
in the past that many different classes of materials 
and systems exhibit memristive behavior including bi- 
nary oxides (Ti02, CuO, NiO, CoO, FesOs, MoO, 
VO2) [2OH25] . perovskite-type oxides (Pri-^Ca^MnOs, 
SrTiOsiCr) [26]^, sulfides (Cu2S,Ag2S) [3T^, semi- 
conductors (Si, GaAs, ZnSe-Ge) [34-36J, spintronics ma- 
terials [37H39J and organics [40-42J. Moreover, hybrid 
memristor-CMOS integrated circuits for reconfigurable 
logic applications as well as memristive memory chips 
combining memristive materials with transistors were re- 
cently developed experimentally [HI [43]. Therefore, fab- 
rication of a memristive processor is possible at the cur- 
rent level of technology. A small-scale version of a mem- 
ristive processor can be built using memristor emulators 

m 

III. NUMERICAL SIMULATIONS 
A. Model 

In this section, we present numerical simulations of the 
dynamics of memristive processors. Numerical modeling 
of memristive networks is easily implemented and thus 
offers by itself a practical computational algorithm for 
maze solving. This computational approach, however, 
requires multiple computational steps (see below) as op- 
posed to the real memristive processor discussed in the 
previous section which needs only a single step to perform 
the whole computation. 

For the sake of clarity, we use a simple model of mem- 
ristor [45] whose memristance (memory resistance) is 
given by 

Rf^ = RoNXij + RoFF (1 - Xij) , (1) 

where Ron and Rqff are minimal and maximal values 
of memristance, Xi^j IS the dimensionless internal state 
variable bound to the region < Xij < 1, and (i,j) are 
grid indexes of a memristor to identify its location in the 
network. Similarly to Ref. i45| we choose the dynamics 
of Xij to be given by 

^=a/.,W, (2) 

where a is a constant and Iij{t) is the current flowing 
through the memristor (ij). At each time step, the po- 
tential at all grid points is found as a solution of Kirch- 
hoff 's current law equations obtained using a sparse ma- 
trix technique. The corresponding change in the mem- 
ristors' states was computed using Eq. (§. 

All numerical results reported in this paper were 
obtained using model parameters Ron = 10 Ohms, 
Rqff = 100 Ohms, Rff{t = 0) = 91 Ohms, x{t = 
0) = 0.1 and a = 10^/(s-74) for all memristors. The 
applied voltage was selected to be equal to 50V during 




FIG. 4: (Color online) Solution of a multiple-path maze. Net- 
work state at t = 0.047 s. The maze solution contains two 
common segments (red dots), and two alternative segments 
of different lengths close to the left bottom corner (blue and 
green dots). The memristance in the shorter segment (blue 
dots) is smaller than that in the longer segment (green dots) 
since the current through the shorter segment is larger and, 
consequently, the change of the memristors' state along this 
segment is larger. The arrow at the bottom indicates where 
a modification of the maze from Fig. [3] has been made (posi- 
tion of the removed segment of the wall) . The resistance is in 
Ohms, the voltage is in Volts and the current is in Amperes. 



the first 0.1 s time interval, and -50V at t > 0.1s. The 
sign of the applied voltage was changed in order to bet- 
ter represent the maze solution as discussed below. Two 
types of mazes were considered: a single-path maze, and 
a multiple- (two-) path maze. 



B. Results 

We have applied the computational scheme described 
above to several mazes. In all cases, the correct maze 
solutions were found. Examples of our calculations are 
shown in Figures |3] and ^ These mazes have been 
mapped on a n x n square memristive processor with 
n = 30. Here, gives the total number of grid points. It 
is evident that the square geometry selected for this par- 
ticular case is not generally required. Mazes of any shape 
(e.g., rectangular, circular, irregular) can be mapped on 
a square memristive processor of sufficient size. 

Figure |3] presents results of solution of a single-path 
maze (see also Supplementary Movie 1). We have found 
that with these parameters it takes approximately 0.06s 
to switch every other memristor along the solution path 
into the low-resistance "ON" state. The resulting se- 
quence of memristors in the low-resistance state repre- 
sents the maze solution as shown in Figure [3^. This maze 
solution can be better seen if we change the sign of the 
applied voltage and wait some time. Then, the resistance 
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of memristors - along the solution path - that are in the 
"ON" state will increase, and the resistance of memris- 
tors in the "OFF" state will decrease. Figure [3]3 captures 
a specific moment of time when these memristances are 
equal. At this moment of time, every memristor along 
the solution path is in the same intermediate state and 
thus the solution of the maze is better visible. 

A two-path maze whose solution is given in Figure |4] 
was obtained from the maze shown in Figure [3] by remov- 
ing a single segment of the wall (its location is shown by 
an arrow in Fig. [4|. In this maze, the current fiow- 
ing through a common segment (red dots in Figure [4| 
splits between two possible paths according to their re- 
sistances. Therefore, memristors in the common segment 
change their resistance faster than those in the two pos- 
sible paths (see Figure [4] and also Supplementary Movie 
2). Moreover, comparing memristors from two possible 
paths, memristors in the shorter path change their resis- 
tance faster than those in the longer path. As discussed 
previously, this allows sorting all possible solutions ac- 
cording to their length in such a way that the resistances 
of the memristors along shorter paths are smaller. 

We note that the type of memristors and simulation 
parameters (such as a and V) have been selected for the 
sole purpose of demonstration of our idea. The resulting 
switching time (of the order of 0.05s, see the above dis- 
cussion) is much longer than the typical switching times 
of nanoscale memristive devices that can be as short as 
5ns (see, for example, parameters of nanoionic resistive 
random access memories in Refs. [12, 46j). Since with 
an appropriate choice of voltage/current magnitude the 
switching time of all the memristors along the solution 
path is on the order of the switching time of a single mem- 
ristor, we can argue that the minimum time required to 
solve a maze of arbitrary complexity by our method can 
be as short as few nanoseconds or even less, since, as 
discussed above, full switching is not required to find a 
solution. Moreover, only one step (a single voltage pulse) 
is needed to find the solution. Therefore, the approach 
we suggest is a priori more efficient than any multi-step 
algorithm of present use. In addition, in simulations with 
memristors described by Eqs. ([T])-([2| we have always ob- 
served the system evolution toward a unique (for a given 
maze) stable solution. Therefore, the success rate of cor- 
rect solution is 100% in our scheme. 

Regarding the time complexity of the numerical mod- 
eling algorithm, it is mainly determined by the solu- 
tion of a set of linear Kirchhoff's current law equa- 
tions. The best theoretical estimate for a linear system 
(the Coppersmith- Winograd algorithm [47J) is 0(n^"^^^). 



This estimate thus provides the theoretical time complex- 
ity of the numerical maze solving algorithm suggested in 
this paper. We note that the time complexity of the 
breadth-first search algorithm that can be used to find 
the shortest path is only slightly better. In particular, 
it is given by 0(V E) ^ O(n^), where V and E is the 
number of vertices and edges, respectively [48j. However, 
we would like to emphasize once more that the hardware 
implementation of the memristive processor requires a 
single computational step, and thus outperforms all ex- 
isting maze solving algorithms. 



IV. CONCLUSION 

In conclusion we have shown how a network of mem- 
ristors - a memristive processor - can solve mazes in a 
massively-parallel way. This approach can be realized 
experimentally with available systems and devices, or 
simply implemented on a computer. The hardware im- 
plementation of the memristive processor is superior to 
any existing maze solving methods and therefore it is 
ideal when the complexity of the maze increases with in- 
creasing local connectivity of the graph. Although we 
have considered a processor based on memristors, a net- 
work of memcapacitors or meminductors [llj can also be 
used for massively parallel calculations. As of now, sev- 
eral experimental systems exhibiting memcapacitive and 
meminductive properties are known [12j. Electronically, 
memcapacitive and meminductive circuits can be emu- 
lated using memristors [49, 50j. 

Moreover, we anticipate that the memristive processor 
can facilitate the solution of - or solve - many other com- 
putational problems. Examples of such problems include 
the traveling salesman problem, graph theory problems, 
etc. The memristive processor can then be used as a com- 
plete computational device, or as a supplemental tool 
for traditional computing hardware. Since memristors 
(as well as memcapacitors and meminductors) are gen- 
erally asymmetric devices, unidirectional graphs can also 
be easily realized on appropriate memristive processors. 
We thus envision their use in a large set of applications 
in both basic science and technology. 
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Supplementary materials to the manuscript Solving 
mazes with memristors: a massively-parallel approach by 
Y. V. Pershin and M. Di Ventra 

Supplementary Movie 1 

This movie shows dynamics of maze solving for a single- 
path maze. The time is in seconds, resistance is in Ohms, 
voltage is in Volts, and current is in Amperes. The 
red line links network grid points that are not separated 
by maze walls and are connected by memristors whose 
Rff < 80 Ohms. Note that at t = 0.1s the polarity of 
applied voltage changes. 

Supplementary Movie 2 



This movie shows dynamics of maze solving for a 
multiple-path maze. The multiple-path maze (see Fig. 
3 of the article) is obtained from the single-path maze 
(Fig. 2 of the article) by removal of a single wall seg- 
ment (see Fig. 3 for details). The time is in seconds, 
resistance is in Ohms, voltage is in Volts, and current is 
in Amperes. The red line links network grid points that 
are not separated by maze walls and are connected by 

memristors whose Rff < 80 Ohms. Note that at t = 0.1s 

''J 

the polarity of applied voltage changes. 



The movies can be found at 
http:/ /www. physics. sc.edu/ pershin/memdevices.htm 



